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Abstract 



■ We review Shuvaev's transformations that relate off- forward parton distributions (OF- 

! PDs) to the so-called effective forward parton distributions (EFPDs). The latter evolve like 

conventional forward partons. We express nonforward amplitudes, depending on OFPDs, 
^3 ' directly in terms of EFPDs and construct a model for the EFPDs, which allows us to consis- 

tently express them in terms of the conventional forward parton distributions and nucleon 
form factors. Our model is self-consistent for arbitrary x, ^, /u, and t. 

1 Introduction 

^ , The treatment of nonforward high-energy processes, such as deeply virtual Compton scattering 

^ I (DVCS) and hard exclusive electroproduction of vector mesons, in perturbative QCD gives rise 

to a new class of parton distributions, the so-called skewed parton distributions (SPDs) |l|, ^, 
^, |, ^, §], 0) §1 and double distributions (DDs) 0, ^, generalizing the well-known conventional 



parton distributions [T^, |TT], and, at the same time, the nucleon form factors. In the following 
we restrict ourselves to the off-forward parton distributions (OFPDs) introduced by Ji 0, ||, |[, 
which are equivalent to nonforward 0, ^ or off-diagonal parton distributions. Therefore, our 
results can be easily generalized to nonforward and off-diagonal parton distributions. 

The off-forward parton distributions Hp{x,C,,t, which parametrize nonforward matrix el- 
ements of light cone bilocal operators (P'|Op(— n/2,n/2)|P)|„2=0) depend on the momentum 
fraction x of the average nucleon momentum P := {P + P')/'^, which the initial state parton p 
carries, on the "skewedness" ^ = —n ■ A/2n ■ P with A := P' — P, on the momentum transfer 
invariant t := A^, and on the renormalization scale fi. For vanishing ^ and t, they are identical 
to the usual forward parton distributions. Detailed reviews on off-forward parton distributions 
can be found, e.g., in Refs. |^. 

Recently, Shuvaev ||13[ demonstrated that the off-forward parton distributions can be related 
(at least in leading order) by simple transformations to so-called effective forward parton distri- 
butions (EFPDs), the renormalization scale dependence of which is governed by the conventional 
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forward evolution equations. These relations have led to some progress in determining the shape 
of the off-forward parton distributions for small values of ^ [Q, since the EFPDs can be identified 
with the usual partons for small values of ^ and arbitrary scale /i. 

In the present paper, we express nonforward amplitudes directly in terms of effective forward 
parton distributions. Furthermore, we define a family of self-consistent models for EFPDs, in 
which the effective forward parton distributions are obtained from the conventional forward parton 
distributions and nucleon form factors at arbitrary scale /i. 

In the next section we briefly review the basic properties of the off-forward parton distribu- 
tions, and we define the effective forward parton distributions. In Sec. ^, we recalculate Shuvaev's 
inverse transformations, which relate the EFPDs to the OFPDs, and we derive their support in 
X, which is not identical to — 1 < x < 1 as for the conventional forward parton distributions]^ In 
Sec. H, Shuvaev's transformation is brought into a form that is convenient for a further analytical 
and numerical treatment. In Sec. ^, we connect the effective forward parton distributions directly 
to nonforward amplitudes, and we briefly discuss the reliability of simple approximative formulas. 
In Sec. 1^, we introduce our model. Finally, in Sec. |^, we summarize our results, and we draw the 
conclusions. 

2 OfF-forward and effective parton distributions 

The long-distance behavior of hard scattering processes, which is not calculable in (QCD) per- 
turbation theory, is factorized in matrix elements of light-cone bilocal operators. A Fourier 
transformation of diagonal matrix elements results in the conventional quark and gluon densi- 
ties q{x) and g{x). Analogously, the off-forward parton distributions are defined by nonforward 
matrix elements: 



{p\s'\i,,{-i)ig i,,{%)\p,s)\ 



v?=0 



= U{P\ S') i U{P, S) J e-^^("-^) H^{x, e, t) dx + 0{A) , (la) 
{P', S' |F;,(-f ) n'^n-^g,, F^\(f )| P, S) 



^2=0 
"+1 



= ^ U{P', S') i U{P, S){n-P) J ^ e-^^("-^) H^{x, e, t) dx + 0{A) , (lb) 

where Q^ab) is the Wilson gauge link. Note that we have an additional factor of x in the definition 
of the gluon distribution compared to the original definition of Ji ^ = x Hg\ which 

removes an "artificial" singularity for finite ^ [§] . Due to time reversal invariance and hermiticity, 
Ji's OFPDs are even functions of ^, so it is sufficient to treat only positive values of ^. The 
different 0{A) contributions can be found, for example, in Refs. ^ |15|. For vanishing A the 
off-forward parton distributions reduce to the diagonal partons: 

i/^(x,0,0)=g(x), (2a) 
Hg{x,0,0)=xg{x) . (2b) 



^We use throughout parton distributions with both signs of x, i.e., q{—x) ~ — q{x) and g{—x) = — g{x). 
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The renormalization of the defining operators leads to a scale dependence of the off-forward 
parton distributions. The evolution of the OFPDs takes a simple form at the one-loop level for 
the Gegenbauer moments, 

GU^,t,ii) := j ^ CCi'/'\l)H^{x,^,t,ft)dx, (3a) 

^niU^^'y■=^^^^^^ £\"-'C^:i?{l)H^{x,t,f,,Odx, (3b) 

since they evolve exactly as the Mellin moments in the diagonal case p[. For example, for the 
Gegenbauer moments of the nonsinglet off-forward quark distributions one has 

/ a (u)\ 

Grit /^) = J /io) , (4) 

where 7on and /Sq are the leading coefficients of the nonsinglet anomalous dimension and the beta 
function. This fact allows the definition of effective forward parton distributions, whose Mellin 



moments equal the corresponding Gegenbauer moments of the OFPDs [13 



J a;"g5 j(x,/i)dx = G^(^,t,/i) , 



(5a) 



j ^ 9^,tix, /i) dx = (^, t, /i) . (5b) 

Their scale dependence is governed by the conventional evolution equations, and they reduce to 
the diagonal quark and gluon densities for — > 0. The crucial point is that the effective and 
off-forward parton distributions can be related to each other by Shuvaev's transformations ||TB|. 
As these transformations do not depend on the momentum transfer invariant t and the scale n, 
we can safely skip them in the following. 

3 Shuvaev's inverse transformation 

We start with connecting the effective forward parton distribution to the off-forward ones. This 



can be done by Shuvaev's inverse integral transformation [|I3 

f+i 



(l^{x) = j !^q^{x,t,y) H^{y,Ody , (6a) 
g^{x) = £\;\x,^;y)H^{y,Ody . (6b) 



We briefly sketch the main steps of the derivation of the integral kernels /C l{x,^; y) in Ref. |T3 



in order to determine the support properties of the EFPDs. The calculation is based on the 
formal inversion of the Mellin moments in Eq. (B|) 



IdiscV^, (7a) 



TT ^-^ X 

n=0 
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Figure 1: Circles of convergence of the power series in Eq. ^ for positive w. The corresponding 
figures for negative w are received by mirroring at the vertical axis. The thick lines show the 
discontinuities of the generating functions. 



X] 



^ oo 

— disc > 



cm 



n=l 



X' 



n+1 



(7b) 



with 



disc F{x) = — lim [F{x + ie) — F{x — ie)] . 

The Gegenbauer moments G^'^(a;) are defined in Eq. (|^). The factorial functions in Eq. 
are replaced using the integral representation of the beta function {Eqs. (6.1.18) and (6.2.1) in 
Ref. [|T6l}. We obtain 



X 



+1 

-1 
+1 



- — disc [ - ^ 



-1 



-— disc 

TT 



l=5:cp.(|)(fi)"d. 

n=0 ^ ' 

n=0 ^ ^ 



(Sa) 



(8b) 



The expressions in square brackets are the integral kernels ]Cgg{x,^; y). Before we state their final 
form, we have to look at the generating functions of the Gegenbauer polynomials {Eq. (22.9.3) 
in Ref. 1T§}: 



exp(-z/Log(l - 2wz + z^)) = J^^n^ 



(9) 



n=0 



The generating functions on the left-hand side analytically continue the power series on the right- 
hand side to the complete complex plane. In Fig. we show the circles of convergence of the power 
series and the discontinuities of the generating functions, which arise from negative arguments of 
the complex logarithm. We see that we have to distinguish two cases: \y\ > ^, which corresponds 
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Figure 2: Support area of the effective forward parton distributions. The dashed lines correspond 
to X = ±'C/2- The contribution of the meson-wave-function-like part of the OFPDs is restricted 
to the region between the dashed lines. 



to the parton-distribution-like region of the OFPDs, and \y\ < ^, which corresponds to the meson- 
wave-function- hke region. Let us begin with the latter case. One might think that for \y\ < C, 
one has no contribution to the integral kernels )C~g {x, ^ ; y) because the generating functions are 
analytical for any finite and real z. However, the discontinuity ai z = oo produces delta functions 
and their derivatives 5^'^\x). From the circle of convergence we find that their contribution is 
in any case restricted to |a;| < ^/2 in the effective forward parton distributions. For |?/| > the 
discontinuity can be easily calculated. We face strong singularities in the generating functions at 
the end points of the cut, therefore, we have to take derivatives of less singular functions, so that 
the s- integral in Eq. (RI) is convergent. An examination of the circle of convergence shows that 



9? I 



X) 



and g^{x) = , if |x| > Xf, — ^ (l + - • (10) 



This defines the support area of the EFPDs, which is shown in Fig. Finally, the complete 
result for the integral kernels of Shuvaev's inverse transformation is 



Ai\ dy 



+ ^(1 - If I) E(-1)'^7^;7^^"C'(3^'^(|) '^^"H^) , (11a) 

n=0 



Ax e{-s^ + s^-^ - f)VT^s 



+ 0(1 - If I) ^^^^^^^ ■ ^^^^^ 

n=l ^ 



We cannot carry out the derivatives, which would lead to infinite contributions from the deriva- 
tives of the theta functions and divergent integrals. The infinite sums represent the contribution 
of the meson-wave-function-like part of the off-forward parton distributions and do not appear 



in Ref. |13[, since the discontinuity at z = oo, resp. x = 0, was overlooked. The occurence of 
delta functions and their derivatives in the effective forward parton distributions are comparable 
to meson-exchange-type contributions in double distributions ^ . 
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The expressions for the integral kernels ]C~g{x,C,', y) show that the Gegenbauer moment inver- 
sion is not practicable, in general. Therefore, it can generally not be used for a simple solution 
of the evolution equations of off-forward parton distributions. 

4 Shuvaev's transformation 

As we will later see, the predictive power of the formalism lies in relating the effective forward 
parton distributions to the off-forward ones by Shuvaev's integral transformation 

Hq{x,0 = j >^q{x,^]y)q^{y)dy , (12a) 

Hg{x,0 = lC,{x,i-y)g^{y)(iy . (12b) 



The full derivation of the integral kernels lCq^g{x,^; y) can be found in Ref. [T^. We merely state 
the finite result: 



^.(^'e;y) = ^#^ r^(^-l)J \. , . d., (13a) 



Ti\/\y\dy J[y\ J-i ^' \ y{l - s"^) ~ X + is 




Again, performing the derivatives would give divergent integrals and infinite contributions from 
the end points. The derivatives of the theta function give rise to the "suspicious overall sign" that 
is mentioned in Ref. . Equation (|I^) is equivalent to previous results presented in Refs. |T^, |TB 



It is useful to express the integral kernels in terms of standard elliptic integrals, because it is 
then possible to perform the derivatives analytically. First, we give the symmetry properties of 
the integral kernels: 

J^q{x,-i;y) = ICq{x,i;y) , ICq{x,i;-y) = +ICq{-x,i;y) , (14a) 

!^g{x,-i;y) = ICg{x,i;y) , ICg{x,i;-y) = -ICg{-x,i;y) . (14b) 

Of course, /Cg,g(x, ^; y) obey the fundamental (— ^ ^ ^)-symmetry of Ji's off-forward parton dis- 
tributions. The relations on the right-hand side show that it is sufficient to restrict the calculation 
to positive values of y. We define 

a := - - ^ - — - Ayx + ^ , &:=-- — + —^4^/2 - 4yx + ^2 (15) 
i 2y 2y^ ^ y > ^ 2y 2y^ ^ y > K J 

which correspond up to the x/^, which we added for convenience, to the zeroes of the denominator 
in the square root in Eq. (|TB|). As the radicand has to be positive, we must further restrict the 
possible range of y to 

y > Xa := ^{x + \/x'^ - i'^) , if x > ^ . (16) 
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With help of the integral tables in Ref. ||T9[ we obtain our final result: 



+ eiy - x,)^, {R^iO, f , 1) - 1^ RniO, f , 1)) for x > ^ , (17a) 



.^(^ + (^HO, ^, 1) - m Rn{0, ^, 1)) for X < e 

X Rp{0, f , 1) - 1'^^$^ Rn{0, f , 1)) for X > e , 

x(f^^HO,^,l)-|^^i?z.(0,^,l)) forx<e, 

where Rp and i?D are Carlson's elliptic integrals of the first and second kind (see, e.g., [[ 

1 r 1 



RF{x,y,z) 



RD{x,y,z) 



2 70 ^/it + x){t + y){t + z) 



dt 



{t + z)^{t + x){t + y){t + z) 



dt 



(17b) 

) with 
(18) 

(19) 



These integral kernels accumulate the main properties of the x and ^ dependence of off-forward 
parton distributions. For x ^ ^, the OFPDs essentially look like forward quarks and gluons: 



1 



]C^{x,^-y) = 6{y-x) + -Oi^,) , 

X 

lCg{x,^;y) = x6{y-x) + 0{^) . 



(20a) 
(20b) 



The forward evolution concentrates the effective forward parton distributions at ?/ ~ 0. There- 
fore, the small-y behavior of the integral kernels }Cq^g{x,^]y) reproduces the well-known asymp- 
totic forms of the off-forward valence and singlet quark, and gluon distributions 



i/:(x,0 = ^(l-|l 



2 /•£ 



H'q{x,O = 0il- 

H(x,i) = e{i~ 



,^^ x{e-x^) 





2\2 /•£ 



- + 0{l)jq^^{y)dy 



(21a) 



^'(f + (^(|)) 1/^7^(1/) . (21b) 



Additionally, these equations prove that the integrals in Eq. (|T2|) are well defined and convergent. 

Even the physical interpretation of the different regions in x holds in the formalism of EFPDs. 
Inserting Eq. (|l^ in Eq. (0) yields 



Hg{x,0= / {}Cq{x,^;y)q^{y) + }Cg{-x,^;y)q^{-y))dy 
Jo 

Hgix,^)= / {}Cg{x,^;y) + ICg{-x,^;y)) g^{y)dy . 
Jo 



(22a) 
(22b) 
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For |a;| > ^ the off-forward parton distributions are related to corresponding effective forward 
(anti)partons with a minimum momentum \xa\- For |x| < ^ the picture of a meson wave function 
is supported by a simultaneous contribution of effective forward partons and antipartons with 
any momentum y. The different expressions for the integral kernels )Cq^g{x,^]y) for |a;| ^ ^ show, 



analogous to |T^, that the off-forward parton distributions are not analytic at |a;| = ^. The 



analyticity of the OFPDs for x ^ ^ requires that the effective forward parton distributions need 
to be analytic for |a;| > only. 



5 Nonforward amplitudes and EFPDs 

The use of off-forward parton distributions is required in deeply virtual Compton scattering 
(DVCS) and hard exclusive electroproduction processes. Detailed information can be found in 
Ref . , and references therein. Here, we are only interested in the part of the amplitudes that 
refers to the OFPDs: 

+1/1 1 



-1 
+1 



MO ■= / r-- + , , . H^i^, dx , (23a 



X — ^ + ie X + ^ — ie 



AiO ■■= / 7-- + • - Hgix, dx , (23b) 



where we have neglected the 0{A) contributions, analogous to Ref. The imaginary part 

of the amplitudes is related to the diagonal elements Hgg{C,,C,), which can be expressed by the 
effective forward parton distributions: 

-nH^^{^,0 = ImAqiO = - / 4VT^ql{^r^^)dz , (24a) 

-—Hg{^,0 = ImMO = -J^ 32z'VT^g^{^^^)dz , (24b) 

with the quark singlet g|(x) = q^{x) — q^{—x). The imaginary part is essentially dominated by 
the behavior of the EFPDs around x ~ ^/2. For small values of C, this region can be accurately 
described by 

xq^ix) ~ x'^" , (25a) 
xg^{x) ~ x"^" . (25b) 



If we insert Eq. (g^) into Eq. (^4|) and set the upper integration limits to one, we achieve approx- 
imation formulas for the imaginary part of the amplitudes in Eq. (^3]): 

(26a) 

g|(|) r(A, + 3) ^ ' 



_ ImA(0 ^ 8v/^r(A, + |) 

9,ii) r(A, + 4) 



R'r ■■= -^^^^ - ■ (26b) 
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Figure 3: Comparison of the exact imaginary part to the approximative ratio R^^g{Xg^g) for q^{x), 



A similar ratio was already presented in Ref. |14|, where, however, the imaginary part was com- 



pared to diagonal partons at x = 2^, which leads to an extra factor 2^+^'^'''9. 

In Fig. 1^ we show a comparison of the exact ratio, derived from Eq. (0), to the approximation 
in Eq. (|26|) for the effective distributions in Eq. (p5[) . The change of the integration limit has 
no remarkable effect up to ^ ~ 0.1. The accuracy of the gluon ratio is slightly worse compared 
to the quark ratio, because of an additional factor of z"^ in the integrand in Eq. {^). Since the 
quotients of the gamma functions in Eq. ( PB] ) have a weak Xq^g dependency and Fig. ^ shows a 
good stability under a change of Xq^g, we can conclude that Eq. (^) is an excellent approximation 
of the imaginary part of the amplitudes for the values of ^, where the effective forward parton 
distributions can be reliably identified with the conventional forward quark and gluon densities. 

The calculation of the real part is straightforward but tedious [we have used Eq. (|13|) rather 
than Eq. (p!7|)]. The principal value integration can be performed exactly and the final result 
consists of integrals without any strong singularities: 



Re Agio = i (-^^ " ^ + 4v/rT^ + aVi^^ g^{^) dz 
+ r 4(i-8z + 4v/^(l + ^))r7^(4)d^ . 



(27b) 



We note that the expressions in brackets in the first integrals in Eq. (|27|) are always negative, 
therefore, the two integrals partly cancel each other. Again, we insert the small-x behavior of 
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(a) Rf%Xg) : i?^^' (exact) 



(b) i?f (Ag) : i?f (exact) 



Figure 4: Comparison of the exact real part to the approximative ratio R^^g{Xg^g) forq^{x), 



X) ~ 



2j -^13,9 1 



Eq. (|25D into these integrals and set the lower limits of the second integrals to zero, so that 
everything can be evaluated and yields the following ratios between the real and imaginary parts: 



R. 



Re 

'<},9 



tan ■ 



vrA 



Q,9 



:28) 



This is identical to the result, achieved by dispersion relations, in Refs. [^, From Fig. |^ 
we see that the quality of the approximations of the real parts is significantly more sensitive — 
note the logarithmic scale — to the change of the integration limit and to a variation of A^ g. 
Additionally, the right-hand side of Eq. (^) depends strongly on Xg g, and the first integrals in 
Eq. ([27| ) have dominant contributions from two regions: around x ~ and x ~ ^/2, i.e., we must 
require that Xg^g is essentially constant for small x. Therefore, only for very small ^, when Eq. ( |25| ) 
is a valid approximation for a large range of x for the usual forward quark and gluon distribution, 
the latter can be used to reliably predict the real part of the amplitude. Nevertheless, for small 
values of A„,„, where the real part is strongly suppressed, the absolute values of the amplitudes 



are determined to a good precision for small ^ [T^ . 



6 Moment- diagonal models 



In this section, we try to build a model for the effective forward parton distributions. The 
situation does not seem to be very promising, because Shuvaev's inverse transformation cannot 
generally be used and we face a difficult support area for the EFPDs in Fig. ||. 

References 521 were dealing with a model for off-forward parton distributions, the 
Gegenbauer moments of which are independent of ^. Because the corresponding EFPDs are 
also independent of ^, these models manifestly violate the support area in Fig. |], as was also 
recognized in Ref. [ff^, though it should be a good approximation for small values of ^. This 
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model would have had the great advantage that it would have been stable against a change of 
the input scale, since the evolution of the Gegenbauer moments is identical to that of the Mellin 
moments. 

Nevertheless, this idea can be used to find valid models for the effective forward parton 
distributions. Because a common n-, or t-dependent factor does not have any influence on the 
evolution equations, we can generalize the model with ^-independent Gegenbauer moments to a 
class of moment-diagonal models with a common proportionality factor: 

Gr (e, t, /i) := const(n, e, t) X (yu) • (29) 



The t dependence is usually factorized, const (n, ^) x Fi{t), where Fi(t) is the Dirac form factor ||Tj, 



These models allow, provided that the Mellin inverse can be performed and leads to valid 
supports, a direct and simple calculation of off-forward parton distributions and amplitudes, with 
help of the formulas (P!?]), (Plj), and (P?]), for arbitrary x, ^, fi, and for the region of t, where 
the 0{A) contributions can be neglected. The simplest model for proton OFPDs of this class, 
which fulfills all known theoretical constraints for t = (e.g., see Ref. f^) and gives a good 
approximation of the t-dependence, is 

(e, t, /i) := 2 (0 ^ T„+i(r^) Ff{t) Mr (/i) , (30) 

where M^'^ (/x) are the Mellin moments of the usual quark and gluon distributions in the proton, 
Ff{t) is the Dirac form factor of the proton with Ff{0) = 1, and Tn{x) are Chebyshev polynomials 
of the first kind {Eq. (22.2.4) in Ref. It is advantageous to use the Gliick-Reya-Vogt 1998 



(GRV 98) parton distributions [jTO[, since they are given for x values down to 10 ^, which allows 
us to compute the real part in Eq. (p7[) for small values of ^ with a high accuracy. With use of 
Eq. (22.3.25) and (4.4.27) in Ref. |T§, a Mellin inversion yields 

?,(^) = - \x\) qi^^^) + 0C-4^ - \x\) , (31a) 

9,{^) = [(^C-^4^-\^\)9{^^^)+0C-^^-\x\)g{j^^ . (31b) 

We see that the effective forward parton distributions are a simple combination of two rescaled 
forward parton densities with the correct support area and an appropriate common t-dependent 
factor. The first one gives the conventional forward quark and gluons for vanishing ^ and t. The 
contribution of the second summand is restricted to < ^/2, i.e., it influences only the meson- 
wave-function-like region |x| < ^ of OFPDs, and is negligible for small values of ^. Therefore, 



most of the numerical results of Refs. ^ |T8| can be accurately transferred to the model in 
Eq. (llD. 

It is an interesting fact that the argument of the first forward parton density in Eq. ( pT] ) is 
very similar to the Georgi-Politzer ^-scaling variable ||23|, ^ = 2x^/(1 + \/l + 4x|M^/(5^), that 



originally described target mass effects in deep inelastic scattering. Hence, one can argue that 
the arguments of the parton densities in Eq. (^) reflect skewedness effects. But such arguments 



can as well be relics of the Gegenbauer polynomials that appear in the derivation p4|, ^ of the 
,^-scaling variable. 
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7 Summary and conclusions 

In this paper, we presented with Eqs. ([T^), P^), and ( P?]) simple expressions that relate effec- 
tive forward parton distributions, which evolve like conventional forward partons, to off-forward 
parton distributions and nonforward amplitudes. We emphasized that the off-forward parton dis- 
tributions and nonforward amplitudes can be directly determined from the conventional forward 
parton distributions and nucleon form factors at arbitrary scale /i for moment-diagonal models. 
Exemplary, we stated a simple self-consistent model for the EFPDs of the proton in terms of 
the GRV 98 parton distributions and the Dirac form factor of the proton, which allows us 
to predict off-forward parton distributions and nonforward amplitudes for arbitrary x, ^, /i, and 
(not to large) t. These predictions should not differ too much from results of other models at 
least at small ^, as the results in Refs. [|T1|, ^ show. 



Nevertheless, it would be illuminating if further self-consistent moment- diagonal models exist, 
especially models that have a qualitatively different behavior in the meson-wave-function-like 
region, such as the off-forward parton distributions of chiral soliton model calculations p6|, since 
the real part of nonforward amplitudes is dominated by this region and gets important for large 

Because of the complicated support area of EFPDs, an investigation of the double distributions 
of moment-diagonal models might be helpful. 
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